arXiv: 1508.06903v2 [physics.plasm-ph] 5 Dec 2015 


Nonlinear Landau damping and modulation of electrostatic waves in a nonextensive 

electron-positron pair plasma 

Debjani Chatterjee^ and A. P. Misra^’j^ 

^Department of Mathematics, Siksha Bhavana, 

Visva-Bharati University, Santiniketan-731 235, West Bengal, India 

Abstract 

The nonlinear theory of amplitude modulation of electrostatic wave envelopes in a collisionless electron-positron (EP) pair 
plasma is studied by using a set of Vlasov-Poisson equations in the context of Tsallis’ q-nonextensive statistics. In particular, 
the previous linear theory of Langmuir oscillations in EP plasmas [Phys. Rev. E 87, 053112 (2013)] is rectified and modified. 
Applying the multiple scale technique (MST), it is shown that the evolution of electrostatic wave envelopes is governed by 
a nonlinear Schrodinger (NLS) equation with a nonlocal nonlinear term oc P J — ^') [where P denotes the 

Cauchy principal value, f is the small-amplitude electrostatic (complex) potential, and ^ and r are the stretched coordinates in 
MST] which appears due to the wave-particle resonance. It is found that a subregion 1/3 < q < 3/5 of superextensivity {q < 1) 
exists where the carrier wave frequency can turn over with the group velocity going to zero and then to negative values. The 
effects of the nonlocal nonlinear term and the nonextensive parameter q are examined on the modulational instability (MI) of 
wave envelopes as well as on the solitary wave solution of the NLS equation. It is found that the modulated wave packet is 
always unstable (nonlinear Landau damping) due to the nonlocal nonlinearity in the NLS equation. Furthermore, the effect of 
the nonlinear Landau damping is to slow down the amplitude of the wave envelope, and the corresponding decay rate can be 
faster the larger is the number of superthermal particles in pair plasmas. 
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I. INTRODUCTION 


Electron-positron (EP) pair plasmas are ubiquitous and play important roles in many astrophysical situations such 
as the early universe [T], Van Allen radiation belts, and near the polar cap of fast rotating neutron stars [5], black 
holes [3], pulsars [3], quasars [5], active galactic nuclei [5], accretion disks |3 as well as in laboratories [5]. In black 
holes, pulsars and quasars such EP plasmas are emitted in the form of ultra-relativistic winds or collimated jets by 
some of their most energetic objects. Because of the intrinsic and complete symmetry between the positively charged 
(e.g., positrons or positive ions) and negatively charged (e.g., electrons or negative ions) particles, the dynamics of 
pair plasmas become significantly different from that of electron-ion plasmas or from a purely electronic beam. Over 
the last few years, a number of experiments have been performed to create EP plasmas (See, e.g.. Refs. 131 and [Tn|l . In 
such experiments, it has been observed that because of the fast annihilation and the formation of positronium atoms, 
the identification of collective modes in EP plasmas is practically impossible. To resolve this issue and to identify the 
collective modes properly, a number of experiments have been proposed to create pair-ion plasmas (See, e.g.. Ref. [TT| 
an dini). However, most recently, ion free high-density (~ 10^® cm ^) neutral EP plasmas with unique characteristics 
have been produced in the laboratory by using a compact laser driven setup [8] . It has been reported that because of 
their unique characteristics together with the charge neutrality, small divergence as well as high average Lorentz factor, 
such EP plasmas can exhibit collective behaviors and thereby opening up the possibility of studying the collective 
dynamics of EP plasmas in a controlled laboratory environment. 

The collective oscillations of EP plasmas and associated wave dynamics together with the formation of solitary 
waves and shocks have been extensively studied over the past two decades (See, e.g.. Refs. IT3HT3|) . In other works, 
the modulational instability (MI) and the nonlinear evolution of electrostatic and electromagnetic wave envelopes have 
also been studied in EP plasmas (See, e.g.. Refs. [T7Ul9)l . However, most of these works are based on hydrodynamic 
models. Using the kinetic theory approach, i.e., using Vlasov-Poisson equations, the effects of Landau damping on 
various kinds of wave modes (linear theory) as well as on electrostatic and electromagnetic solitary waves [governed by 
the Korteweg-de Vries (KdV) equation] have also been studied in EP plasmas [H] and in some other environments (See, 
e.g. Refs. [2nH23l and references therein). However, to our knowledge, no theory for the explanation of electrostatic 
wave modulation and the effects of nonlinear Landau damping on MI and nonlinear evolution have been reported in 
nonextensive EP plasmas. 

Many spacecraft observations (e.g., in the Earth’s bow-shock, upper ionosphere of Mars, the vicinity of the Moon 
etc.) |24j and laboratory experiments |25j confirm the presence of nonthermal and superthermal particles which do not 
follow the Maxwellian distribution, but show some deviation from the thermodynamic equilibrium. The presence of 
such nonthermal particles has also been confirmed in astrophysical environments |2()) . Several models for phase-space 
plasma distributions with nonthermal or superthermal wings or other deviations from Maxwellian distribution have 
been proposed in recent years. One of such distributions, which was first reported by Renyi m and subsequently 
proposed by Tsallis [35], is the Boltzmann-Gibbs-Shannon (BGS) entropy in which the degree of nonextensivity of 
the plasma particles is characterized by the entropic index q. The distribution function with q < 1 characterizes 
the system with more superthermal particles (superextensivity), while the distribution with q > 1 indicates plasmas 
containing a large number low-speed particles (subextensivity) compared to the Maxwellian one {q 1). Such q- 
nonextensive distribution has been widely considered in a number of works to investigate various collective plasma 
modes and nonlinear coherent structures (See, e.g., Refs. [T31133H3T1I . It is to be noted that the k distribution function 
and the g-nonextensive distribution function in Tsallis’ statistics are somewhat equivalent in the sense that in both 
these cases, the spectrum of the distribution functions show similar behaviors. In fact, there is a formal transformation 
K = 1/{1 — q) which can provide the missing links between these two velocity distribution functions |32] 

On the other hand, it is well known that waves in plasmas can undergo collisionless damping when they resonantly 
interact with trapped and/or free particles, i.e., when the particle’s velocity approaches the wave phase velocity 
or group velocity. Such collisionless damping was first theoretically predicted by Landau [33], and later confirmed 
experimentally by Malmberg et al. [33j. Motivated by these inventions, Ott et al. [3D] first theoretically investigated 
the effects of linear Landau damping on the nonlinear propagation of ion-acoustic solitons in electron-ion plasmas 
through the description of KdV equations and on the assumption that particle’s trapping time is much longer than that 
of Landau damping. Later, Ikezi et al. [35j . based on their experimental investigation, have emphasized the importance 
of nonlinear Landau damping in plasmas. Accordingly, Ichikawa |3bj investigated the the effects of nonlinear Landau 
damping due to resonant particles having the group velocity of the wave, on the modulation of electrostatic wave 
envelopes in electron-ion plasmas. In this work, he assumed that the typical time scale is much longer than the 
bouncing period of particles trapped in the potential trough. It has been shown that the nonlinear wave-particle 
resonance leads to the modification of the nonlinear Schrodinger (NLS) equation with a nonlocal nonlinear term 
proportional to a Gauchy principal value integral as well as the local (cubic) nonlinear term. Eurthermore, in Ref. 
1361 it was reported that, in contrast to the ordinary NLS equation, the nonlinear resonance always leads to MI of 
wave envelopes against a plane wave perturbation regardless of the sign of PQ positive or negative, where P is the 
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coefficient of the group velocity dispersion and Q is the cubic nonlinear term of the NLS equation. 

In this work, our aim is to revisit and extend the work of Ichikawa |36j in an unmagnetized collisionless EP pair 
plasma in the context of Tsallis’ g-nonextensive statistics. Starting from a set of Vlasov-Poisson equations and using 
the reductive perturbation technique (RPT), we show that the NLS equation is not only modified by the nonlinear 
resonant effects but also by the nonextensive parameter q which contributes to both the dispersive and nonlinear 
(local and nonlocal) terms. It is found that the nonextensivity significantly modifies the wave frequency, the group 
velocity, the nonlinear frequency shift and the energy transfer rate for the modulated wave packets as well as the 
amplitude of solitary wave solutions of the NLS equation. 


II. BASIC EQUATIONS AND DERIVATION OF NLS EQUATION 

This section mainly focuses on the derivation of the NLS equation. Though, the relevant analysis is almost the 
same as in Ref. 1361 we, however, review the analysis and present the subsequent derivations and explanations for 
the expressions of various plasma modes in more details. We consider the nonlinear propagation of electrostatic wave 
packets along the x-axis in a collisionless electron-positron pair plasma. The basic equations for the dynamics of EP 
plasmas are given by the following Vlasov and Poisson equations 

dFg dF^ _ e^d^dF^ ^ 
dt ^ dx rua dx dv 


^ J F^dv, (2) 

where Fa is the distribution function with its unperturbed part Fa\v), v is the particle velocity and (/> is the 
electrostatic potential. Also, the particle charge and mass are given by = —e, TOq = nie for electrons and = e, 
rua = rUp for positrons. We assume that the equilibrium state (at t = 0) is spatially uniform field-free EP plasma and 
that the perturbation from the equilibrium state is purely electrostatic. Furthermore, we consider the equilibrium 
distribution of electrons and positrons f1^\v) to be the ^-distribution function as in the Tsallis’ nonextensive statistics 
P5] . It is to be mentioned that due to the resonance of plasma particles having the group velocity of the wave, the 
distribution functions become singular, and so the direct application of RPT to the Vlasov-Poisson equations Q 
and ([^ rnay not determine uniquely the contributions of the resonant particles. In order to treat these singularities 
properly we introduce the multiple space-time scales as |36j 

X ^ X + e~^ri + t^t + e~^a. (3) 

We assume that the amplitude of the carrier wave (with wave number k and the wave frequency ui) is infinitesimally 
small, and so for t > 0 a slight deviation o(e)] from the uniform equilibrium value will occur. Thus, we expand 


OO OO 


Fa{v,x,t) = F^°\v) ^ fj^^{v,ri,aX)e^p[il{kx-ujt)] 


n—1 1— — 00 


OO OO 


(4) 


= Y (’7,o-,Oexp[i/(A;a:-a;f)], 


TL — l 1 — — 00 


where e is a small (< 1) positive parameter measuring the weakness of the wave amplitude and 

4>^-i = for reality conditions to hold. Here, the asterisk denotes the complex conjugate quantity. We note 

that the expansion @ for Fa is valid only in the non-resonance region where |u — u}/k\ ^ o(e) is satisfied. In the 
resonance region where v « w/fc, the above expansion may not be appropriate to apply to the Vlasov equations. So, 

(n) (n) 

the components ( and (j)] ’ are further expanded into multi-scale Fourier-Laplace integrals as 

-I n poo 

= -77rF2 dKf^^^{v,K,n,C)exp[i{Kr)-na)] 

(^TT)^ Jc J-ao ’ 

1 p poo 

</>|"^(»7,cr,C) = / dfl dK^['^\K,nX)exp[i{Kr]-na)], 

Jc J-oo 

where the contour C is parallel to the real axis and lies above the coordinate of convergence. 
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We first derive the evolution equation of electrostatic wave envelopes for arbitrary species of particles with unper¬ 
turbed arbitrary velocity distribution function, and then we extend our analysis in neutral EP plasmas in the context 
of Tsallis’ 5 -nonextensive statistics [5H]. So, we substitute the stretched coordinates @ and the expansions @ into 
Eqs. 0 and 0 to obtain, respectively, as 




rr?_ ^ ^ ^ ^ I I jji 


S = 1 l' — — 00 


oo oo 

E E 

l' = — oo 


_9 


(n-s-1) 


d 


71 — 8—2 

h-v 


JLas) 


( 6 ) 


- 2ilk-^y - i2lk^y - 47r^ e„ J f^dv = 0. (7) 

where we have used the symbol = to denote the equality in the weak sense and disregarded the terms which contain 
and in Eq. Q. 

In the subsequent analysis, we determine the contributions of the resonant particles having the wave group velocity 

(n) (n) 

by solving the cr-evolution of the components ( and (j)i ' as an initial value problem with the initial condition 

/a”oE)^>o'= 0>C) = 0, n>l, (8) 

in the multiple space-time scheme corresponding to that on the distribution function 

/a’}](^,^ = 0)=0. (9) 


A. Harmonic modes for n = I — 1: Linear dispersion law 

Equating the coefficients of e from Eqs. 0 and Q for n = 1, / = 1, we obtain (For details see Appendix [A| the 
following linear dispersion law: 


fc -I- 47r "y^ Ca 

together with the linear Landau damping rate given by 


Go,{v) 
LO — kv 


IL = 




dv = 0, 


Go 


Ic {uj — kv) 


.dv, 


( 10 ) 


( 11 ) 


where denotes the Cauchy Principal value, and we have made the analytical continuation of the integral over v 
along the real axis passing infinitesimally above and under the pole at v = uj/k with the constraint of weakly damped 
waves and 


Ga{v) 


6a d 
ma dv 




( 12 ) 


B. Harmonic modes with 1^0, n = 1: Some conditions 


From Eqs. ^ and Q, we equate the components for I ^ 0 and n = 1, and then use the dispersion relation (10) to 
obtain the following conditions (See Appendix [B| 


y\ = 0 and = 0 for |^| > 2. 


(13) 
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C. Zeroth harmonic modes for n = 1, 2; / = 0 


Here, we examine the second order terms with n = 2 and I = 0. Thus, we have from Eq. 


— -I- 1,— _ C — 0 


while the first order terms with n = 1 and Z = 0 of Eq. Q yields 

fi]ldv = 0- 


( 14 ) 


(15) 


Substituting the Fourier-Laplace integrals given by Eq. ([^ into Eq. (14) and solving it as an initial value problem 
we obtain 


K, H, C) = - K,a = 0,C)- H, C), 

where we have used an arbitrary constant as —f^Q {v, K,a = 0, (). 


(16) 


Next, a substitution of Eq. (16) into Eq. (15) gives 




{K,fix) = 


= 0,C) 

n — Kv 


dv/'^t 


G„ 


n — Kv 


dv. 


(17) 


Using the initial condition (|^ and taking Fourier inversion of Eqs. (16) and (17) we obtain for n = 1, 1 = 0 the 
following zeroth order components: 


f(l) ^ Q 

Ja,o ~ 


a(1) 


= 0 . 


(18) 


The vanishing of and hence of in fact eliminates the contributions of the modes associated with the singu¬ 
larities of the integrals in Eq. 0- Thus, in this way the initial condition ([^ uniquely defines the present problem 
to investigate nonlinear auto modulation of (w, k) modes. 


D. Modes with n = 2, 1 = 1: Compatibility condition 

We proceed to examine the the second order, first harmonic modes with n = 2 and I = 1. Thus, from Eqs. (§ and 
([^ we obtain the following compatibility condition (See for details Appendix [C|) 

+ = (19) 


where A is given by 


A = 


4^5]( 


Ga 


Ic (w - kv)'^ 


dv 


1 -|- 47r 5^ I 


Ic (w - kv)'^ 


G„dv 


( 20 ) 


We find that this expression of A is exactly the same as the group velocity Vg = dw/dk obtained by differentiating 
the dispersion relation (10) with respect to k. Equation (19) shows that the a — rj variation of can be related to 
a new coordinate defined by 


such that 


^ = rj — Xa = e(x — At), 


( 21 ) 


( 22 ) 


This indicates that the coordinate ^ in Eq. (21) establishes a clear relationship between the reductive perturbation 
theory and the multiple space-time expansion method. 
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E. Second harmonic modes with n = I — 2 


For the second order quantities with n = I = 2, we obtain from Eqs. Q and Q the following expressions (See 
Appendix 


„( 2 ) ^ _ 

Ja,2 — 


to — kv + iv 


'^a(P2 o™ 


Ga 


2ma dv \uj — kv + iv 


( 1 ) 


(23) 




where 


Aik, w) = 47r 


^ mo 


1 d f G. 


I(y u) — kv dv \uj — kv 


dv. 


(24) 


(25) 


Note that in Eq. ( p^ the contributions of the resonance terms at the phase velocity oj/k are neglected because our 
basic assumption [Eq. 0 ] for EP plasma oscillations is that the linear Landau damping is a higher-order effect than 
second order. 


F. Modes with n = 3, 1 = 0 


Erom Eqs. §) and 0 we consider the terms corresponding to n = 3, 1 = 0, and use the relations ( [T^ and ( |l^ to 
obtain a set of reduced equations which, after use of the Eourier-Laplace transforms with respect to rj and a and the 
initial condition ®, yield 


f (2) 

Ja,0 


= k^ 


W{K,n) K 
A(^){K,n)n- Kv 


nia n — Kv 



H{K,n), 


(26) 



,.2 H{K,n) 

AM (AT, n) 


yv{K,n), 


(27) 


where the relevant details and the expressions for H{K,Q), 'W{K,^l), A*^'^^(Ar, $7) and Ia{v) are given in Appendix 
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G. Third order harmonic modes with n = 3, 1 = 1: The NLS equation 


Considering the terms for n = 3 and 1 = 1, we obtain from Eqs. ([^ and Q a set of reduced equations, which after 
few steps (See Appendix]^, result into the following modified nonlinear Schrodinger equation for the small but finite 
amplitude perturbation 4>{^,t) = as 


'dr^^de 


+ Q\4’\'^4’ 4 —^ 

TT 


c-e 


-f iS(j) = 0. 


(28) 


The coefficients of the dispersion (group velocity), cubic nonlinear (local, nonlocal nonlinear terms, respectively, are 
P, Q and R, given by P = {lj2)d'^ujjdk'^ = P/a, Q = 7 /a and R — S/a, where 


dv, 

O' ^ 

(29) 

Q, V / 

(30) 

/ 1 1 \ n 

(31) 

"'“(eT + p)*’ 
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6 = —^{k,oj) 


( 32 ) 


with 


B{k, Lo) = 47r 





r 1 d 

\ " ^ 

( \1 

ml j 

u! — kv dv 

uj — kv dv 

\uj — kv J 


dv. 


(33) 


Furthermore, the coefficient S of Eq. (281, representing the linear Landau damping rate associated with the resonant 


particles having the phase velocity of the carrier wave, is given by 

g ^ SjphL 


(34) 


In Eq. (28), the coefficient P appears due to the group velocity dispersion of the wave envelope. The most significant 


contribution of the resonant particles having the wave group velocity is the appearance of the nonlocal nonlinear term 
oc R. This resonance contribution also modifies the local nonlinear coefficient Q, which appears due to carrier wave 
self-interactions originating from the zeroth harmonic modes (or slow modes). 


III. CONSERVATION LAWS 


Before we proceed to the modulation of electrostatic wave envelopes in g-nonextensive EP plasmas, we verify some 
important conservation laws that are associated with the ordinary NLS equation (Here, ordinary means in absence of 
any nonlocal nonlinearity in the NLS equation). We will show that the nonlocal nonlinear term in Eq. (28) violates 
these conservation laws. We note that in absence of the Landau damping effects, i.e., for R = S = 0, the NLS equation 
( |28[ ) possesses an infinite number of conservation laws. The first three conserving quantities are namely, the mass 
Ii = f \4>\'^d^, the momentum I 2 = (2i)“^ / and the energy /s = / — {Q/‘^P)\4>\'^) d^- The 

similar quantities, however, for the NLS equation (28) with Landau damping satisfy the following equations: 




/i = 0, 


(35) 


dr 


+ 2 S]l 2 + -P 




(36) 


dr 


25)73+A// 


c-e 




d_ 




didi’ = 0 , 


(37) 


where the terms oc S and R (the Cauchy principal value integrals) are due to the linear and nonlinear Landau damping 
effects. Next, in Eq. ( [T7| using the fact that the integral over is a convolution of the functions V [1/(C^ — 0] and 
where and noting that the Fourier inverse transform of 

i sgn(s) = —{l/Tr)V{l/^), we have 


^ 2 


exp(is^')|s|(^(s, T)ds. 


So, performing the integral over as a Fourier transform of |(/'(^^'^)P we obtain 

V 


d^di' = ^ [ |s|(^(s,r)|0(-s,T)pds, 






(38) 


(39) 


where ‘hat’ denotes the Fourier transform with respect to ^ or . Furthermore, using (p(s,T) = —2is|(()(s, r)p we 
obtain from Eq. (37) 


d 


R 


—-f2S')/3 = -— / s |()>(s,r)| |(/)(-s,t)| ds. 


(40) 


From this equation we observe that if the linear Landau damping is a higher-order effect than , the term cx S can be 
neglected, and so the the left-hand side of Eq. (40) represents the rate of change of the wave energy. Also, the integral 


7 




















on the right-hand side is positive dehnite, and to be shown later that for long wavelength EP plasma oscillations, R 
is always positive for 1/3 < g < 1 and g > 1, and negative for q < 1/3. Furthermore, it has been shown in Ref. 
M and will be shown in this work also that the relevant results may not be valid for q < 1/3. Thus, in both the 
superextensive and subextensive EP plasmas we have the inequality (the equality holds for </> = 0 V 


dh 

dr 


< 0 , 


(41) 


implying that an initial perturbation (e.g., in the form a soliton) will decay to zero with time r, and hence a steady 
state solution with j/aj < oo of the NLS equation (28) may not exist in presence of the nonlinear Landau damping 
term in nonextensive EP plasmas. 


IV. ELECTROSTATIC ENVELOPES WITH NONEXTENSIVE STATIONARY STATE 


The NLS equation (28) in Sec. IIG has been derived in a general way to describe the evolution of electrostatic waves 


in plasmas with arbitrary species a and with arbitrary velocity distribution function for equilibrium plasma state. 
However, the main purpose of the present work is to investigate the dispersion properties of carrier wave modes, the 
linear Landau damping rate associated with the resonant particles with the phase velocity, the modulational instability 
and nonlinear Landau damping due to resonant particles having the group velocity as well as nonlinear evolution of 
wave envelopes in EP plasmas with g-nonextensive stationary states. In the latter, the g-distribution function in one 
space dimension is given by 




1 - (go, - 1) 


rrictV 

2kBTa 


l/(9c-l) 


(42) 


where TOq and Tq, are respectively the mass and temperature of a-species particles (a = e, p, respectively, stand for 
electrons and positrons) and is the Boltzmann constant. The spectral index g^ gives a measure that determines the 
slope of the energy spectrum of nonthermal particles. It also measures the deviation from the thermal equilibrium state 
(go —>■ 1). We consider a fully symmetric and charge-neutral EP plasma in which T^, = Tp = T and rUe = rup = m. 
Such assumptions are justified with the experimental works in pure pair plasmas with particles having the same 
dynamics mm- Also, in the creation of a pure EP plasma, the whole system reaches a common thermal state with 
equal particle temperature. So, the spectral index for electrons and positrons are taken as the same, i.e., qe = qp = q 
M- The normalization constant Aa^q is given by = Lqy/mo,/2TrkBTa = Lq/y^vt, where Vt = y/ksTjm is the 
particle’s thermal velocity and Lq is defined by 


Lq= < 


1 , 9 —>■ 1 , 

Hi-g 2 ) 




I. 


(43) 


Note that in the extensive limit g 
Boltzmann distribution 


I, the distribution function (42) reduces to that of the standard Maxwell- 
(v) = i/rTq/T^TrfcB^exp (—77iQ,u^/2/cBTa). However, in the present work the case of 
g —>■ I may not be recovered directly because of simplihcations of some expressions in the subsequent analysis, e.g., 
for long-wavelength perturbations. The distribution function with g < I represents the superextensivity, i.e., presence 
of more (compared to the Maxwellian) particles with velocities larger than their thermal velocities (superthermal 
particles), whereas the case of subextensive distribution is represented by g > I, meaning that there is a large number 
of particles with velocities lower than their thermal velocities. Furthermore, because of long-range nature of Coulomb 
forces between the plasma particles and the presence of many superthermal particles in astrophysical environments 
[26], a g-distribution with g < 1 is strongly suggested in real plasma systems or superthermal plasmas. For more 
about the behaviors of the distribution function F^\v) for different values of g, the readers are referred to, e.g.. Ref. 


M 

In what follows, we consider the nonlinear propagation of high-frequency (w > ujp) and long wavelength [k <C 

where kd = (STrnQe^/kBT'j is the Debye wave number and n^o = ripo = uq is the equilibrium plasma number 
density] oscillations in g-nonextensive plasmas whose phase velocity greatly exceeds the thermal velocities of electrons 
and positrons, i.e., uJr/k ^ Vt > v. In this case, the linear dispersion laws and the coefficients of the modified NLS 


















equation (281 can be simplified. To this end, we evaluate the Cauchy Principal value integral in Eq. (A4) with the 


following expansion for the integrand 


dv 




dv 


k 

UJr 


(uJr/k) — V 

I dFi°\v) dF^°\v) 


dv 


dv 


' h h'^ t-3 

^ r\j nj <-) rv o 

1 H- V H - ttV^ H -■ 

, UJr UJr Wr 


dv, 


(44) 


where the integration limits are taken as ±Umax = ioo for q < I and zLy/2kBT/m{q — 1) for q > 1. These limits are 
considered due to the fact that for g > 1, the distribution function (421 has a thermal cutoff, which limits the velocity 
of particles to a maximum value, i.e., v < Umax- However, such cutoff is absent when g < 1. In this case, the velocity 
of particles remains unbounded (For details, see, e.g.. Ref. IT^ . 

Next, to evaluate the integrals in Eq. (44) we note that the g-distribution function F^\v) is an even function of 
v, while dFa\v) / dv is an odd function of v. Thus, one can evaluate the integrals as follows: 


Umax dP^^^ (v\ 

v'^—^^^dv = 0, form = 0,2,4, 


dv 


Umax ^^^ 0 ) 

V -- dv = — 1, 


dv 

,dFi°\v) 

dv 


:,dF^^\v)^ 

V -;r- dv = 


dv 


'3g- 1’ 


(45) 


dv = — 


60uC 


(3g-l)(5g-3)- 


We have evaluated the above integrals by parts and the average value of as 




0,,2 

v'^F^\v)dv = ‘ 


3g- 1 


(46) 


We mention that for g > 1, the above integrals, in which the limits are ±Umax, are obtained by reducing the integrals 
to beta functions of the form i?(m,n) with m, n > 0 and finally to gamma functions using the relation B{m,n) = 
T{m)T{n)/T{m + n). Similar integrals can also be evaluated for —1 < g < 1 in which the limits are ±oo by using the 
relation B{m,n) = dx and the above relation between beta and gamma functions. However, in each of 


these cases of superextensive and subextensive plasmas, we will obtain the same results by means of Eq. (43) except 
the factors — g for — 1 < g < 1 and ^g — 1 for g > 1. From Eq. (46) we also note that the values of 1) are 
further restricted to the region 1/3 < g < 1, because otherwise, the average value of may diverge. In particular, 
in the limit g —>■ 1, Eq. ( [4^ reduces to the well known energy equipartition relation {^mv'^) = ^ksT. Thus, our 
results may be valid for both the superextensive (1/3 < g < 1) and subextensive (g > 1) distributions of electrons 
and positrons in plasmas. 


Now, we use t he re sults as in Eqs. (45) and (46), and in the region of small wave numbers (1 >> = fc^/A:^) to 

obtain from Eq. (A4) the following dispersion relation for electrostatic carrier waves in a nonextensive EP plasma: 


CUr - 


1 + 3x2 


+ 60x^ 


3g-l ' (3g-l)(5g-3)J ’ 

together with the linear Landau damping rate, obtained from Eq. ( [n| ), as 


7L = 


TT iUpLq 
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1-(9-1)U-7 


1 


2x2 3g- 1/ 


{2-q)/{q-l) 


(47) 


(48) 


where ujp = is the plasma oscillation frequency in a charge neutral EP plasma. In the same way, the 

group velocity expression (20) reduces to 


A = 


GujpX 


(3g - l)kd 


45g-ll 2 ' 
(3g-l)(5g-3)^ . 


(49) 


The dispersion relation (47) exactly agrees with that obtained in Ref. [Uup to the second order correction term oc x^. 
However, we observe some interesting features by retaining the term oc x^ in Eq. (47) which were overlooked in Ref. 
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M These interesting features, as can be seen from Fig. are that the carrier wave dispersion curve and hence the 
group velocity of the wave envelope can turn over through the y-axis, going to zero values and then to negative values. 
These reduction of the wave frequency and the group velocity occurs in the superextensive sub-region 0.47 < q < 3/5. 
Beyond this region, i.e., for q > 3/5, both tOr and A increase with increasing values of %. Furthermore, it is observed 
that as one goes from the superextensive sub-regime 3/5<g<lto subextensive one with q > 1, the wave frequency 
and hence the group velocity are seen to get reduced. This is expected since more the superthermal particles the larger 
is the phase velocity, in agreement with the previous results M- However, some disagreements are also there for the 
linear Landau damping rate 7 ^. Here, we mention that though its analytic expression is the same as obtained in Ref. 
m however, the features we observe here are quite distinct. In fact, the possibility of growing instability, as shown in 
Ref. El cannot be made in the superextensive or subextensive plasma regions. From the mathematical expression of 
7 l one can check that it is always negative regardless of the values of q and y in 1/3 < g < 1 (or g > 1) and o < y < 1 
respectively. These are clear from Fig. Physically, since the phase velocity of the carrier wave is assumed to be 
larger than the particle’s velocity, the wave modes may lose energy to the particles instead of gaining energy from 
them, and thus be damped. From Fig. we also find that there are two subregions of y, namely 0 < y < yo and 
Xo ^ X ^ 1- In ths former, the damping rate increases, while in the latter the same decreases with increasing values 
of y (See the left panel). Furthermore, the damping rate is seen to be higher in plasmas with more superthermal 
particles and with wavelength in (0 ;< y < yo). However, the same can be true with higher the number of low-speed 
particles (or with increasing values of g) in the regime yo ^ X ^ 1- These are also clear from the right panel of Fig. 

Here, as g increases, the damping rate increases, however, it gets reduced at long-wavelength perturbations (See 
the dashed line). 


(a) (b) 




FIG. 1. (Color online) The real part of the carrier wave frequency Urjuip [Subplot (a); Eq. ( |47| l] and the group velocity X/vt 
[Subplot (b); Eq. ( |49[ )] of the wave envelope are plotted against the carrier wave number y = k/kd for different values of the 
nonextensive parameters as in the figure. It is seen that the carrier wave dispersion curve can turn over with the group velocity 
of the wave envelope going to zero and then to negative values (See the solid and dashed lines) in the subregion 0.47 < g < 3/5 
of the superextensive region 1/3 < g < 1. In the region 1/3 < 1 < 0.47, the values of X/vt greater than unity are inadmissible, 
while in the other regime g > 3/5, both the frequency and the group velocity assume only the positive values. 


Next, we calculate the various terms (in the limit y^ <C 1) which appear in the coefficients of the NTS equation 


(281 for the g distribution (42) as follows: 


2k 

a = -I 1 

UJp 


X 


3g- 1 


X 


*d(3g - 1) 


50g-ll 2 
(3g-l)(5g-3)^ . 


H = 0, W = 0, U = 0, 


B = -15 


(-) 

\m/ 


2 


1-b 


24 


3g- 1 


X 


(50) 

(51) 

(52) 

(53) 


10 































(a) (b) 



FIG . 2 . (Color online) The imaginary part of the carrier wave frequency [The linear Landau damping rate given by Eq. 

( | 48 | )] is plotted against x = k/kd [left panel (a)] and the nonextensive parameter q [right panel (b)] for different values of q 
and X as in the figure. From the left panel it is seen that in one subregion 0 < x Xo, the damping rate increases with x, 
while in the other Xo ^ X ^ 1 h reduces with x. Furthermore, lower the percentage of superthermal particles, the higher is the 
Landau damping rate. The right panel shows that the damping rate becomes weaker in the limit of long-wavelength (x ^ 1) 
oscillations. 


A = - 


2Lgk^ 


2q-3 
81 / 71(1 - q) V 9 - 1 


1 - 


54(3g2-59 + 3)(g-l) 2 


( 3 g-l) 2 ( 2 g- 3 ) 


X 


(54) 




X 


90g2 - 1899 + 97 2 
( 39 - 1 ) 2 ( 59 - 3 )^. 


(55) 


C = 


1 2Lq 49-5 
KbT J 1/71(1 — 9) 89 — 3 


/ 2 ( 2169 ^ - 1779 + 53 ) V 7 r(l- 9 ) 45 ( 9 -l) \ 2 

i (39-1)2(49-5) Lg 49-5 P 


(56) 


D = - 




\KbT J k 


1 /7r\i/2 6 L 0 


(i) 


39- r 


909^ + 219 - 61 2 

( 39 - 1 ) 2 ( 59 - 3 )^, 


Thus, the coefficient P of the NLS equation (|28|) reduces to 

P = 




11 


( 39 -l)fcU ( 39 - 1 )( 59 - 3 )' 


(57) 


(58) 


To reduce the other coefficients, namely Q and R, we first obtain the expressions for 0(fc,w) and $(fc,a;) in the limit 
X^ 1. So, we calculate, respectively, the resonant and non-resonant contributions to 0(/c,a;) as 




(59) 


A ) \kBTj ^^(1-9) 3 ( 1 -9)^ ■ 


2kLg 49 — 5 2 

X 


(60) 


Thus, it turns out that the resonant contribution to the nonlinear coupling coefficient Q being smaller can be disre¬ 
garded to obtain 


Q = \ 


2L„ 


49-5 2 

WpX ■ 


2 \kBTJ 1/7r(l — 9) 3(1 — 9) 


(61) 
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However, the contribution from the group velocity resonance through ^{k,uj) gives rise the coefficient R as 


R = 




3g — 1 yksT 




1/2 




( 62 ) 


Note that the above expressions for P, Q and R are obtained for 1/3 < g < 1. The similar expressions for g > 1 can 
also be obtained by replacing the factor (1 — q) by {q — 1) under the square root in the expressions for A, C and Q. 
In Sec. [Vj we will find that though the condition for the modulational instability does not depend on the sign of PQ, 
but on the presence of R, however, the sign of PQ are important for determining the values of the frequency shift 
and the energy transfer rate (in particular, their maximum values) as well as their values in some particular cases, 
namely when the wave intensity exceeds or smaller than a critical value which depends on P, Q and the wave number 


of modulation. On the other hand, in Sec. Ill we have seen that whether a steady state solution of the NTS equation 


exists or not depends on the coefficient R. Furthermore, in Sec. VI we will also examine the effect of R on solitary 
wave solution of the NTS equation. Thus, it is useful to investigate the properties of P, Q and R. 

Inspecting on the expressions for P, Q and P, which explicitly depend on the nonextensive parameter q and the 


< y < 1), we find that for superextensive plasmas with 1/3 < q 1, we have 
for subextensive plasmas, we have (5>0((3<0)inl<q<5/4(g>5/4) 


nondimensional wave number y (0 
Q < 0. However, when q > \, i.e. 
and for 0 < y < 1. Also, P > 0 in both the superextensive and subextensive regimes, i.e., 1/3 < g < 1 and g > 1 
with 0 ;< y 1. Furthermore, P, which has singularity at g = 3/5, is positive for g > 3/5, i.e., for 3/5 < g 1 

(superextensive) and g > 1 (subextensive). However, in the superextensive subregion 1/3 < g < 3/5, P can be 
negative or positive depending on the rage of values of y in 0 < y < 1 (See the left panel of Fig. [^. The regions for 
P and PQ are shown as contour plots in the y — g planes in Fig. We hnd that in the superextensive regime, PQ 
is negative for 3/5 ^ g ^ 1, however, it can be positive or negative in the sub-regime 1/3 < g < 3/5 depending on 
the values of y in 0 y < 1 (See the middle panel). From right panel (third from left) we find that in subextensive 
plasmas, PQ < 0 for g > 5/4 and for 0 < y < 1 except in some small regions of g (> 5/4) and y in which PQ > 0. 
The latter is also true for 1 < g 5/4 and 0 < y ;< 1. We also note that P, PQ are undefined at g = 3/5. 



FIG. 3. (Color online) The regions for P > 0, PQ > 0 (shaded or gray) and P < 0, PQ < 0 (blank or white) are shown 
in the y — g (y = k/kd) plane. From the left panel (a), it is seen that P > 0 for q > 3/5, i.e., in the regimes 3/5 < g < 1 
(superextensive) and g > 1 (subextensive). However, in the other superextensive regime 1/3 < g < 3/5, it can be positive 
or negative depending on the values of y in 0 < y < 0.1. Middle panel (b) shows that in the superextensive regime, PQ is 
negative for 3/5 < g < 1, however, it can be positive or negative in the sub-regime 1/3 < g< 3/5 depending on the values of y 
in 0 < y < 1. From right panel (c) we find that in subextensive plasmas, PQ < 0 for g > 5/4 and for 0 < y < 1 except in some 
small regions of g (> 5/4) and y in which PQ > 0. The latter is also true for 1 < g < 5/4 and 0 < y < 1. Note that P, PQ 
are undefined at g = 3/5. 


V. THE NONLINEAR LANDAU DAMPING AND MODULATIONAL INSTABILITY 

Here, we follow the same analysis as in Ref. 1361 Though, the relevant analysis is standard, however, we repeat here 
for the sake of readers. It is well known that when the group velocity dispersive coefficient (P) and the local nonlinear 
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(cubic) term (Q) of an ordinary NLS equation have the same sign, i.e., PQ > 0, its plane wave solution exhibits 
instability against a plane wave perturbation of its amplitude and phase [40) . However, the present modified NLS 
equation (28) contains, in addition to the usual dispersive and nonlinear terms, the nonlocal nonlinear term which is 


associated with the resonant particles having the group velocity of the wave. Thus, it is of interest to examine how 
the nonlocal term contributes to the modulational instability and describes the nonlinear Landau damping process. 
Here, we assume that, before modulation, the NLS equation (28) has a plane wave solution of the form 


exp 



where p and tr are real functions of ^ and r. Substituting the solution (63) into Eq. (28) we get 

d 


dr 


9 , , 

P+ 


(63) 


(64) 


9 9 2PR^d 

-a + a-p = 2PQ-p+—P- 


A^de + p^- 




(65) 


As before, we assume that the linear Landau damping term is higher order than and linearize Eqs. (64) and (65) 
by splitting up p and a into their equilibrium (with suffix 0) and perturbation (with suffix 1) parts as: 


P = Po + Pi cos - Qt) + p 2 sin - Hr), 


( 66 ) 


cr = tJi COS {K^ — fir) + (72 sin {K^ — fir), 


(67) 


where H and K are, respectively, the wave frequency and the wave number of modulation. Now, under the perturba¬ 
tions (1^ and (1^, the solution (1^ can be expressed as 


A- cos (kox-ujot) 

-l-Ai(Ar) cos [{ko + eK)x — (wq -I- eXK + e^H)t] 
+A 2 {K) sin [(fco + eK)x — (wq -I- eXK + e^H)t] 
-|-similar terms with K —K and H —>■ —H, 


where 


MK) = 


4 ^ ^ APR 




J = 1, 2. 


( 68 ) 


(69) 


It follows that due to the linearization of Eqs. (64) and (65), the electric potential (j) [Eq. ( [6^ ] describes the three- 
wave interaction of the unperturbed carrier wave (fco,WoTand two side bands with wave numbers and frequencies 
fcoiejATj and wq ± eAjArj ± Now, substituting the perturbation expansions ([^ and (67) into Eqs. (64) and (65), 


and assuming the smallness of the Landau damping coefficient, i.e.. S' = 0, we obtain after eliminating ai and (72 the 
following equation 


-P 2poPQK^ - P^K* 
2poPR sgn{K)K^ 


-2poPR sgn(K)K^ 

+ 2poPQK^ - P^K^ 


Pi 

P2 


= 0 , 


(70) 


where the coupling between pi and p 2 appears due to the nonlocal nonlinear term (the second term on the right-hand 
side) in Eq. (65). Thus, for nonzero values of pi and p 2 we obtain the following dispersion relation for electrostatic 
wave envelopes in collisionless nonextensive EP plasmas. 


(H^ -h 2poPQK^ - P^KY = -{2poPRK^Y. 


(71) 


The negative sign on the right-hand side of Eq. (0 shows that whatever be the signs of P and Q, the Langmuir 
wave packet is always unstable due to the nonzero coefficient R associated with the resonant particles having the wave 
group velocity. Since Eq. 0 is, in general, complex in H, we seek a general solution of it by considering H -PiL 
with Hr, r being reals, and obtain 


Hr = i—^ 

V2 


|(p2^2 _ 2poPqY + {2poPRfy^^ + {P^K^ - 2poPQ) 


1/2 


\Kl 


(72) 


13 


























|(p2^2 _ 2poPQf + {2poPRfy^^ - {P^K^ - 2poPQ) 


-I 1/2 


1 ^ 1 , 


(73) 


where we consider the upper (lower) sign for isT > 0 {K < 0). 

In what follows, we consider two different limits of the wave amplitude. In the small amplitude limit with po ‘C 
\P/2Q\K'^, Eqs. (72) and (|7^ reduce to 




±PK\ V^tpoR. 


(74) 


Since a linear dispersion relation can be expanded as uj{k) = a;(fco)+w'(fco)(A: —fco) + (I/2)w"(fco)(fc —^o)^H-> we take 


the upper sign for if > 0 and lower sign for if < 0 in Eq. (74). Thus, for if > 0, we have 12^. = PK^, E = —poR in 


which the imaginary part (oc R) is solely due to the resonant particles having the group velocity of the wave envelope. 
It follows that the Langmuir wave packets under the modulation is unstable. Since i? > 0 for both superextensive 
(1/3 < g < 1) and subextensive {q > 1) plasmas, this instability is a kind of decay, and of course independent of P 
and Q. Thus, in the small amplitude limit as above, i.e., when the wave intensity po is well below a critical value, 
the real part of fl relates the group velocity dispersion, while the imaginary part of fl describes the nonlinear Landau 
damping process. In the latter, the wave energy is transferred from the higher frequency side bands to lower frequency 
ones. From Eq. (74) we also find that the frequency shift can be positive or negative depending on the values of q 
and X 9'S in Fig. For example, 11^. > 0 in the regimes 3/5 < g < 1 (superextensive) and g > 1 (subextensive). 
However, in the other superextensive regime 1/3 < g < 3/5, it can be positive or negative depending on the values 
of X in 0 < X 0.1. Furthermore, for a fixed wave number of modulation if, the frequency shift gets significantly 
reduced in the region 0.47 ^ g 3/5. This is a consequence of the results to the fact that in this region, the carrier 
wave frequency turn over with the group velocity going to zero values and then to negative values (See Fig. [^. In 
the other region of g, i.e., for g > 3/5, the frequency \ujr increases with x- 

In the large amplitude limit with po |P/2Q|if^, the frequency shift and the transfer rate T can be obtained 
as 




±^/po{-PQ)K^ 


yi+wof 


-I 1/2 


(75) 


r = TVpo{-PQ)K^ 


^l + iR/Qf-l 


1/2 


(76) 


which require PQ < 0. Thus, it turns out that when the wave intensity po greatly exceeds a critical value, the frequency 
shift and the energy transfer rate can be obtained only in the regions of PQ < 0 as in the middle and right panels 
of Fig. both in superextensive and subextensive plasmas. Furthermore, we note that and T are proportional 
to y/oo instead of po as in the small amplitude case. In particular, for i? = 0 or when the cubic nonlinearity (local) 
greatly dominates over the local nonlinear term, the modulated wave becomes unstable for K < = -y/2po|(3/P| as 

in the ordinary NLS equation. 

Next, from Eq. (73) we also find that for a given value of po, the maximum value of the growth rate T can be 
achieved at the wave number Km and for PQ > 0 where 


2_ /Q^ + R^ 

Km -Po( 


The corresponding maximum values of and F are thus obtained from Eqs. (72) and (73) as 


nm = ±^{Q^ + R^f%o, 


(77) 


(78) 


rm = T{Q^ + R^y^%0- (79) 

It follows that the maximum values of the frequency shift and the energy transfer rate for modulated waves can be 
achieved only in the regions of g and x where PQ > 0. Figure]^ (Middle and right panels) confirms that a wide range 
of values of x as well as g, both for superthermal and subextensive plasmas, exist for which PQ > 0. From the above 
results we also conclude that in contrast to the ordinary NLS equation (as in fluid theory) in which the modulational 
instability occurs only for PQ > 0, the modulated Langmuir wave packets in q-nonextensive plasmas always becomes 
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unstable by the effects of resonant particles having the group velocity of the wave irrespective of the sign of PQ > 0 
or PQ < 0. In the former, the maximum values of the frequency shift and the growth rate are achieved for arbitrary 
amplitude of the pump (unperturbed) wave which may not be obtained in Maxwellian plasmas |36j . whereas the latter 
gives asymptotic values of the same for larger values of the wave intensity. 

In a general manner, we numerically examine the properties of fir and F given by Eqs. (72) and (73) for different 
values of the nonextensive parameter q. The results are displayed in Fig. |^for both superextensive and subextensive 
plasmas. From the upper panel of this figure, we find that corresponding to the superextensive regime 0.47 ^ q ^ 3/5 
where the carrier wave frequency turn over with the group velocity (Fig. [^, two subregions of x exist, in one of which 
rir decreases having cutoffs at lower %, while it increases with increasing values of q and x- In the other region of q, 
namely, q > 3/5, the frequency shift is seen to increase with x without any cutoff, while it decreases with increasing 
values of q. For the damping rate F, some different features are observed (See the lower panel). Here, as q increases, 
the values of |F| decrease, i.e., more the number of superthermal particles the higher is the rate of energy transfer 
from high-frequency side bands to low-frequency ones. 



1 


FIG. 4. (Color online) The nondimensional freqnency shift Qr/i^p [upper panel (a)] and the energy transfer rate F /ujp [lower 
panel (b)], given by Eqs. (72) and (731, are shown with respect to the nondimensional carrier wave nnmber y = k/kd in 
superextensive (0.47 < g < 1) and subextensive {q > 1) plasmas. In the snperextensive subregion (0.47 < q < 3/5) where the 
group velocity of the wave envelope vanishes (Fig. [^, the freqnency shift and the energy transfer rates are also seen to have 
cutoffs. 


VI. NONLINEAR LANDAU DAMPING OF SOLITARY WAVE SOLUTION 


In absence of both the linear and nonlinear Landau damping effects, i.e., for R = S = 0, the modified NFS equation 


(28) reduces to the following ordinary NFS equation 

d<j) 




(80) 


In the case where the modulational instability occurs for PQ > 0, a stationary solution (bright soliton) of Fq. (80) 
can be obtained as m 


(/ = i/o exp (i6>), 

where the amplitude (j)q and the phase 9 are given by 

- vqt 


— y 9o 


sech 


1 

TP 


Voi ■ 


L 




(81) 


(82) 


(83) 
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with L(j)Q = y/2\P/Q\ being a constant and z = {^ — vqt)/L denoting the transformation in the moving (with velocity 
Vq) frame of reference. 

In this section our aim is to determine the effect of a small amount of the nonlinear Landau damping (R) associated 
with the resonant particles having the group velocity of the wave envelope on the soliton solution (81). Here, we 


disregard the small effect of the linear Landau damping rate oc S which has been assumed to be higher order than 
e^. From Eq. (41) it is evident that an initial perturbation of the form (81) will decay to zero. Thus, one might 


expect that the amplitude (/)o is no longer a constant but can decrease slowly with time, i.e., (/)o = (f’oiZyT). We 
consider |P|, |(5| 3> |P| ~ e S' ^ p > 0 and do the perturbation analysis of the NLS equation ( p^ with R as 
the small parameter. It can be easily verified that this assumption is valid for both superextensive and subextensive 
regions of q as mentioned before. We follow the similar approach as has been applied in different studies, however, to 
Korteweg-de Vries (KdV) equations (See, e.g.. Refs. andand references therein). Now, under the transformation 
z = (^ — vot)/L, Eq. (p^ reduces to 


del) ■'>Jodcl> P d'^cj) 2 


4,+ Lt, = 

IT J Z — Z' 


(84) 


where we rewrite i? = r to denote R as small (We will replace again r by i? in the final solution) and d(l)ldz' = d(j)/dz 
at z = z'. 


In what follows, to investigate the solution of Eq. (84) we generalize the multiple time scale analysis with respect 
to r, i.e., we consider the solution as 

(/)(z, t) = + • • • (85) 

where i = 0,1,2,3,..., are functions of r = tq,ti,T 2 , ■■■ Substituting Eq. (85) into Eq. (84) and equating the 
coefficients of the zeroth and first orders of r, we successively obtain 




where Ai and A 2 are given by 


A P d -Vo , f d 1 


|<^(o)(zOP 


dz' . 


( 86 ) 


(87) 


( 88 ) 


It can easily be shown that under the boundary conditions, namely /dz, jdz^ —>■ 0 as z —>■ 00 , Eq. 

(86) possesses a solution of the form sech z exp (f 01 z) = 0 some real values of 0i. Now, 

for the existence of solution of Eq. ( p7| , it is necessary that be orthogonal to all solution g{z) of L'^[g\ = 0 

which satisfy p(±oo) = 0, where L'^ is the operator adjoint to L defined by 


/ /)i{z)L[\p2{z)\dz = / 't/i2{z)L'^[f/i{z)\dz, 

J —00 J —00 

with ?/;i(±oo) = ^/> 2 (±oo) = 0, and the only solution of L'^[g\ = 0 is g{z) = sech z. Thus, we have 

/ OO 

A 2 <^^°^sech zdz = 0, 

-OO 

which gives 


J_ df/o 

dr. 


\f^d\ sinh (7r0i/2) 
1 — cosh (ttOi) 


- 

7 T 


coo poo 


( sech^z' 


— OO —OO 


-00 W- 


sech^zexp(i0iz)(izdz' = 0. 


(89) 


(90) 


(91) 


Equation (91) is a first-order differential equation for the wave amplitude c/q. So, its solution can be obtained as 

-1 


00 = 000 1 — i — 

To 


(92) 
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where 0o = <?!>oo at r = 0 and tq is given by (rewriting now r as R) 


Tn ^ = 




00 




cosh {tt9i) — 1 


sinh {ttOi/2) 


^oo ^oo /oppn2 /\ 

V / ( -^ I sech^i; exp(i0i2)(izdz'. (93) 

J— oo J —oo \ ^ ^ J 


Thus, when PQ > 0, an approximate solitary wave solution of the NTS equation (28) with a small effect of the 
nonlinear Landau damping is given by 


= \/</>oo sech z exp(i0), 


where 6 = + (fl — v^ji) t] /2P. 


On the other hand, when PQ < 0, a stationary solution (dark soliton) of Eq. (80) can be taken as [4T 


(j) = (I)oexp{i0), 


where (po and 9 are different from those given by Eqs. (82) and (83), i.e.. 


4>o = (po tanh 


^-VqT 


(94) 


(95) 


(96) 


^ 2P 


vo^ + hpQpl - y I "T 


(97) 


with Lpo = y/\2P/Q\ = a constant. 

Proceeding in the same way as above up to Eq. ( |87[ ), it can be easily verified that tanh z exp {i92z), for 

92, is a solution of Eq. (86) if and only if dp^^^Jdr = 0. Also, for the existence of a solution of Eq. (87) we have 


some 


the same necessary condition but with different g{z) = tanh z. The condition gives 


/ OO 

tanhzdz = 0, 

-OO 


(98) 


from which we obtain 


i\/2tt —— 

OTl 


m 


02 sinh (7r02/2) 
1 — cosh(7r02) 


-pnV 


/ tanh-^ z' 
/-oo V ^ 


tanh^ z exp(i92z)dzdz = 0, 


(99) 


where po = poo at r = 0. As before, the solution of Eq. (99) in the case of PQ < 0 is given by 

- 1/2 


Po — Poo 1 ~ * 


To 


where (Rewriting r as R) 


Tn-^ = I - 


3/2 


Rplo [1 “ cosh(7r6<2)] 


i5(r)(l — cosh7102) + O 2 sinh (^) 


1'°° f tanh^ z' 

/-oo \ Z- z' 


tanh^ z e:xjp{i 02 z)dzdz , 


( 100 ) 


( 101 ) 


with S{t) denoting the Dirac delta function. Thus, when PQ < 0 an approximate solitary wave solution of the NLS 
equation (|28[) with a small effect of the nonlinear Landau damping is given by 


' = 9^>oo 1 ~ * 


To 


- 1/2 


tanhz exp(/0). 


( 102 ) 


where 9 is given by Eq. (97). 


with a smal 


Erom Eqs. (94) and (102), it is evident that the absolute value (|(/)|) of the wave amplitude decays slowly with time 
effect of the nonlinear Landau damping. Figure exhibits a qualitative plot of the absolute value of 


p given by Eq. (94) in plasmas with superextensive and subextensive velocity distributions. It is seen that faster the 
decay of the wave amplitude, the larger is the percentage of superthermal particles in plasmas. 
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FIG. 5. (Color online) The nondimensional decay rate \(i>o/4>oo\ is shown for both superextensive (solid and dashed lines) and 
subextensive (dotted line) plasmas. It is seen that faster the decay of the wave amplitude, the larger is the percentage of 
superthermal particles in plasmas. 


VII. CONCLUSION 


We have investigated the amplitude modulation and the nonlinear evolution of electrostatic wave envelopes in a 
collisionless electron-positron pair plasma in the context of Tsallis’ q nonextensive statistics. Starting from a set of 
Vlasov-Poisson equations and applying the reductive perturbation technique, the dynamics of the wave envelopes 
is shown to be governed by a NLS equation with a nonlocal nonlinear term arising from resonant particles having 
the group velocity of the wave envelope. Such wave-particle resonance also modifies the local nonlinear (cubic) 
coupling coefficient of the NLS equation. Furthermore, the nonextensive parameter q, which measures the excess of 
superthermal particles in plasmas, is shown to modify the dispersive (group velocity), local nonlinear (cubic) as well 
as the nonlocal nonlinear terms significantly. An interesting effect of q is that a subregion (0.47 ^ ^ 3/5) of the 

superextensivity (1/3 < q < 1) exists where the carrier wave dispersion curve {uJr) can turn over with the group 
velocity (A) going to zero and then to negative values (Fig. [^. Such features of the dispersion curve have not been 
reported in the previous work m where the same plasma model has been considered to investigate linear Landau 
damping of Langmuir oscillations in q-nonextensive EP plasmas. In Ref. [T4]these effects were absent due to truncation 
of the wave frequency up to Furthermore, our results in the linear theory show that the electrostatic wave is 
always damped due to resonant particles having phase velocity of the wave in both superextensive and subextensive 
regimes and also for long-wavelength perturbations. These are also in disagreement with the results of Ref. [HI where 
the possibility of growing instability has been predicted. We, however, stress that such a growing instability should 
not appear for high-frequency oscillations in g-nonextensive EP pair plasmas. This is due to the fact that since the 
wave phase velocity is assumed to be larger than the particle’s velocity, the wave modes can be damped by losing their 
energy to the particles. It is found that for a fixed q, two subregions of y exist, in one of which the linear damping 
rate 7 l becomes higher, while in the other it gets reduced with increasing values of y. On the other hand, for a fixed 
y, the absolute value of 7 ^ increases with increasing values of q, however, its value gets reduced in larger wavelength 
of perturbations. 


In the nonlinear regime, we have verified the conservation laws as applicable for a NLS equation. It is found that 
unlike the ordinary NLS equation, the nonlocal nonlinear term, associated with the nonlinear wave-particle resonance, 
violates the conservation laws, leading to a decay of the wave amplitude with time and thereby forbidding the existence 
of a steady state solution of the modified NLS equation (28). We show that the modulated wave packet is always 
unstable (regardless of the sign of P and Q) due to the nonlocal nonlinear term which can describe the nonlinear 
Landau damping process in which the wave energy is transferred from higher to lower frequency sidebands. 

The frequency shift (11^) and the energy transfer rate (F) for the modulated waves are also examined by the 
parameter 5 in a general way as well as in the limits of small (po ^ Pc) and large (po Pc) amplitudes, where 
Pc = \P/2Q\K'^ is some critical value of the pump wave intensity po. It is found that both fir and F attain their 
maximum values only when PQ > 0. Such maxima of Qr and F may not exist in pair plasmas or electron-ion plasmas 
with Maxwellian distributions [36) . However, they assume some asymptotic values in the limit po S> Pc and when 
PQ < 0, which are oc instead of po as in the opposite limit po <C Pc- The regions of q and y for which P, PQ > 0 
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and P, PQ < 0 are also obtained in both superextensive {q < 1) and subextensive {q > 1) regimes. 

The general expressions of both fir and T are also studied by the effects of q. It is found that corresponding to 
the superextensive regime 0.47 ^ q 3/5 where the carrier wave frequency vanishes with the group velocity (Fig. 
[^, two subregions of x exist, in one of which a significant reduction of fir occurs having cutoffs at lower y, while 
it gets enhanced with increasing values of q and y. The existence of such cutoffs are the consequences of the turn 
over effects of the carrier wave frequency as well as the group velocity of the wave envelope. In the other region of q, 
namely, q > 3/5, the frequency shift is seen to increase with y (having no cutoff), while it decreases with increasing 
values of q. However, quite distinct features are observed for the energy transfer rate T. Here, as q (> 1/3) increases, 
the values of |r| decrease, i.e., more the percentage of superthermal particles the higher is the rate of energy transfer 
from high-frequency side bands to low-frequency ones. 

We have also studied the effects of a small amount of the nonlinear Landau damping (oc R) on the solitary wave 
solution of the ordinary NLS equation by assuming that the linear Landau damping oc xl rate is higher order than 
e^. It is found that the wave amplitude decays with time and the decay rate can be faster the larger is the number 
of superthermal particles in EP plasmas. To conclude, the results should be useful for the evolution of nonlinear 
wave envelopes and associated wave damping in collisionless pure EP plasmas or pure pair-ion plasmas such those in 
laboratory [5^, space Pll and astrophysical environments. 
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Appendix A: Modes with n = Z = 1 


Equating the coefficients of e from Eqs. (§ and (0 for n = 1, Z = 1 we successively obtain 


For nonzero perturbations, Eq. (All yields the following linear dispersion law: 

'■ G^{v) 


fc -f 47r I 


' — kv 


dv = 0. 


(Al) 


(A2) 


The small parameter e, which measures the amplitude of the wave, can be related to the linear Landau damping rate 
7 l of the plasma wave as xl ~ , where p is a nonnegative integer and ujr is the real part oi lo = Wr + i^L- Next, 

from Eq. (10) we have the plasma dielectric function 


D{k, w) = 1 — 


47 re^ f dv 


F^°\v)+F^^’iv) 


( 0 )/ 


mk"^ 


V — uj/k 


^ = 0 , 


(A3) 


where for pair plasmas or EP plasmas We = nip = m. If xl AC lOr, i-e., 7 l ~ then the dielectric function 

D{k,uj) can be Taylor expanded in the smallness of xl to obtain its real (with suffix r) and imaginary (with suffix i) 
parts as m 


Dr{k, LOr) = 1 — 


n ih \ /47re2\ 

Di(k,Wr) = -TT I 1 


47re^ 

r dv 

Ff\v)FF'i'\v) 

J 

Q V — LUrjk 


d_ 

dv 


{Fi°)(u)+FW(^;)} 


(A4) 

(A5) 


U—LOrjk 


Next, neglecting the terms of order and higher in Eqs. (A4) and (A5), the expressions for LUr and xl can be 

obtained from the following relations 


7L = 


Dr{k, LUr) = 0, 
Di{k,LUr) _ TT ■ 


dDr{k,LUr)/duJr 






Ic {lu — kv) 


rdv. 


(A6) 

(A7) 
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Appendix B: Modes with n = 1, 1^0 


We equate the components for Z ^ 0 and n = 1. Thus, from Eq. ([^ we have 

il{uj — kv)f^\ + ilkGa(t>\^'^ = 0- 

This gives 

fa] = - 

LJ — kv + lU 


(Bl) 


(B2) 


where v = ±|i^| for I ^ 0 has been introduced to anticipate that the solution in the linear approximation decays with 
Landau damping rate. 

Similarly, for Z ^ 0 and n = 1 equating the coefficients of e from Eq. Q we obtain 


fk'^(j)\^^ J faj dv = 0. 


Next, substituting Eq. (B2) into Eq. (B3) and noting that 


f cZu ~ f dv±iTT f dvS{u! — kv) ^ ( dv±o{e^'^^), 
J R J C 'J 'J C 


we obtain 


Pk + 47r 




(B3) 


(B4) 


(B5) 


where the contour C enables us to remove in the denominator of Eq. (B5). Thus, by means of the dispersion 
relation (10), one must have 


fa] = 0 = 0 1^1 — 2 - 


(B6) 


Appendix C: Expressions for n = 2 and I = 1 

We consider the second order, first harmonic modes with n = 2 and Z = 1. Thus, from Eq. ([^ we have 


i{uj - kv)f^^\ + ikCaCji^i^ = ^fa} + ^§^fa} “ 


which gives 


A2) ^ _ 
Ja,l — ■ 




uj — kv + iv 

Also, from Eq. ([^ we have for n = 2 and Z = 1 the following equation 

- 2ik-^(j)^^^ - 47ry^ Cq, J fa}dv = 0. 


A -k (‘2) 

This equation can be reduced, by substituting from Eq. (C2) the expression of ( into it and using the dispersion 
relation (A2| and the expression (B2), to 


(Cl) 


(C2) 


(C3) 


JL/i) 

drj^ 




Ga 


Ic (w - fcu)2 


d ,(1) d (1) 


dv = 0, 


Equation (C4) can be written in the form 


d . d \ ( 1 ) 


-+A-|<^i (,?,a;C) = 0. 


(C4) 


(C5) 
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Appendix D: Second order harmonic modes for n = I — 2 

For n = I = 2 we have from Eqs. ([^ 

2i{uj - kv)f^l + 2ikGa4>2^ = 




dv"" 


which by means of Eqs. (B2) and (B 6 ) yields 

A2) ^ k 

uj-kv + iv 

d 2 ) ; 




Next, the expression for (^2 obtained from Eq. 0 for n = / = 2 and using the relation (|D2[) as 


i.e., 


- 47 r^ Ca J faldv = 0 , 

(t>2^ = , 


(Dl) 


(D2) 


(D3) 


(D4) 


Appendix E: Harmonic modes with n = 3, 1 = 0 


For n = 3 and I = 0 we have from Eq. 


A/-(2) , A2) iLA2) _ _jfc_e5LA2)_^ Ai) , 


Too 


au+ 977^-1^7;^“4+ m„ 977^1 ^ ’ 


where we have used the relations (18). This equation can further be reduced after few steps using Eqs. (B2), (C2) 


and (C5) to the following equation: 

la + -g^/ri - 


^ "(2) , 9 ( 2 ) ^ d ( 2 ) ^ ea 2r ^ 1 ^( 1 ) |2 


dr]' 


dr] 


(E2) 


where 


d 


= o:: 


'-A 


dv [ (w — fc 7;)2 


Gr 


(E3) 


Next, taking the Fourier-Laplace transform of Eq. (E2) with respect to 77 and a and using the initial condition 
we obtain 


(A-.n,o - e^^^UvmK.n). 


where H{K,n) is defined as 


(j)[^\r]-Xa,C) l^= / dKH{K,n)exp[i{Kr]-na)], 


with 


H{K,n) = 2TTS{n- KX) J + K'). 

Now, for n = 2, 1 = 0 Eq. ^ reduces to 

dTT^Ca f = 0 . 


(E4) 


(E5) 


(E 6 ) 


(E7) 
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Thus, we obtain from Eq. (E4| the following relation 


,2 H(K,n) 


= k 

0 AW(i4:,n) 


where 




w{K,n), 


K 


= -4nKj2 


ITT'Ot J {OJ — — Kv) 

f G, 


Gadv, 


n — Kv 


dv, 


Next, a substitution of Eq. (E8) back into Eq. (E4) results a slow beat wave solution given by 

K 


f(2) ^ p 

Ja,0 — 


W{K,n) K e, 

-Ga{v) - 


A^<^){K,n)n- Kv 


VL — Kv 


Ia(v) 


H{K,n). 


(E8) 

(E9) 

(ElO) 

(Ell) 


Appendix F: Modes with n = 3, 1 = 1 and the NLS equation 

We can now use the expressions for the lower-order quantities so obtained to determine the components for n = 3 
and / = 1. Thus, we have from Eq. (rol) 




a(3) 


UJ 


-kv^^ 


, Gg 9 , 

UJ — kv dr] 


<,( 2 ) 


+ i 


di) 
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(w — kv) 
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G^ d 
UJ — kv d^'^^ 
d d 


Ja,l ^ Q^Ja,l 


da 


‘ dt] 


iTia UJ — kv dv \uj — kv'*'^ 2ma UJ — kv dv \uj — kv) v^i J v" i rria uj — kv 
This expression for is then substituted in Eq. 0 for n = 3, Z = 1 to obtain the following equation 
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k^ d 
:v c 
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-f AdTT V ^ 
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J {uj — kv)^ \ da dt] 

^ f i d _ 

1Tla J UJ — kv dv \uj — kv 

el, f 1 9(1 d 
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E + fid + ^ A fid _ G E^(i) t dv 




G„ 


9?7" 


'9ry' 


UJ — kv dv \ UJ — kv dv \uj — kv 
1 




^ p(2) 


UJ — kv dv 




dvkrjj^i'i = 0 , 


Jy UJ — kv 


(FI) 


(F2) 


where 9(p) is unity for p = 0 and vanishes otherwise, and T denotes the path of integration around v = uj/k va 
the anticlockwise sense. Note that this integral results from the term fa idv — fi^ldv representing the Landau 
damping term. 


Now, substituting the expression (Ell) into the coefficient of of the sixth term of Eq. (F2) we obtain the term 
as 

E (j / jdfa; =i‘‘ I I 

cx ^ ' 

X \w{K, nf/Aid {K, ft) + C(K, fl)l , (F3) 
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where 


C(K.n)=4,J2^ 


K 


Ia{v, K, Vt)dv. 


7c ~ kvY{Vl — kv) 

By deforming the contour C to the Landau contour with real Q we have 

1 1 ,1/17' 


— kv + iiy Q, — kv \K\ \ K 


(F4) 


(F5) 


Thus, using the relation (E6) for H{K,fl), the functions defined by Eqs. (E9), (ElO) and (F4) can, respectively, be 
expressed as follows: 


with 


Similarly, 


with 


W{K,n) = -W{k,uj;X)-tUik,u;;X)^ 

\K\ 


VF(fc,a;;A) =47ry 

f ^ m _ 


1 Gc 


Wq, 7 (w — kv)"^ V — X 
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Also, we have 


where 


D{k,uj-X) = An . 


V- J„(A) 


TO^ (w — kXy 
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AW(i^,f}) = A(A) + *r(A)^ 


A(A) = 47r y Ca f dv, 


(Fll) 
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(F13) 


r(A) =47r2y ecG„(A). 

a 

We note that the group velocity A = duj/dk is independent of K and ft and we have the relation 

^JdnJ dKe.p^iKr,-na)^^H{K,n)=^'^ j 


(27 


in terms of the constrained coordinate (211. Thus, we finally obtain the term given by Eq. (F3) as 

w ,(i) 


MttV ^ 
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m„ / UJ 
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where the symbols used are 


&{k,uj) = 


A r 

-(Vt/2 _ jj2^ 2- WU + C 


(F17) 


$(fc,a;) = k^ 


A2 + r2 


(^2 _JJ2 _2 


A 


A2 + r2 


WU + D 


(F18) 


Here, the functions U, D and F represent contributions of the resonant particles having the group velocity of the wave 
envelope in EP plasmas. Now, the third term in Eq. (F2) can be reduced by eliminating with the help of Eq. 


(B2| as 


-dTT^, 


v-X „ , ( 1 ) 


(u,- 


{u) — kvY ^ dr]‘^ ^ ’ 


(F19) 


in which the first term of it cancels the second term of Eq. (F21 by means of the expression of the group velocity (201. 


Furthermore, inspecting on the first term of Eq. (F2) we find that the derivative Xd/d^ can be transformed into the 
derivative d/dr by redefining the variable Q = Xt. Finally, substitutions of Eqs. (24) and (|F16 ) into Eq. (F2) give 
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